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^ ■ Abstract 

We present a general formalism for extracting information on the fundamen- 
tal parameters associated with neutrino masses and mixings from two or more 
. ^ ' long baseline neutrino oscillation experiments. This formalism is then applied 

^ ■ to the current most likely experiments using neutrino beams from the Japan 

Hadron Facility (JHF) and Fermilab's NuMI beamline. Different combina- 
tions of muon neutrino or muon anti- neutrino running are considered. To 
extract the type of neutrino mass hierarchy we make use of the matter effect. 
Contrary to naive expectation, we find that both beams using neutrinos is 
more suitable for determining the hierarchy provided that the neutrino en- 
ergy divided by baseline (E/L) for NuMI is smaller than or equal to that of 
JHF. Whereas to determine the small mixing angle, #13, and the CP or T 
violating phase 8, one neutrino and the other anti-neutrino is most suitable. 
We make extensive use of bi-probability diagrams for both understanding and 
extracting the physics involved in such comparisons. 
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I. INTRODUCTION 



The solar neutrino puzzle, which has lasted nearly 40 years, has finally been resolved 
by KamLAND [1] after invaluable contributions by numerous solar neutrino experiments, 
especially recently by SNO [2] and Super-Kamiokande [3]. The Large-Mixing- Angle (LMA) 
region of the Mikheev-Smirnov-Wolfenstein (MSW) [4] triangle [5] has now been uniquely 
selected. This resolution to the solar neutrino puzzle has opened the door to the experimen- 
tal exploration of CP violation in the Lepton Sector. Together with the existing evidence 
for neutrino oscillation that has been obtained by the pioneering atmospheric neutrino ob- 
servations [6] and the long-baseline accelerator experiments [7], it lends further support to 
the standard three-flavor mixing scheme of neutrinos. 

The remaining task toward the goal of uncovering the complete structure of the lepton 
flavor mixing would be a determination of the (1-3) sector of lepton mixing matrix which 
is now called as the Maki-Nakagawa-Sakata matrix (MNS) [8]. A full determination of the 
(1-3) sector includes a determination of #i 3 , the CP or T violating phase 5 and the sign of 
Amy. (We define Am^ = m 2 - — m 2 , where rrii is the mass of the i-th eigenstate.) 

The sign of Am 2 3 signals which pattern of neutrino masses nature has chosen, the normal 
hierarchy, Amf 3 > 0, or the inverted hierarchy, Amf 3 < 0. This hierarchy must certainly 
carry interesting information on the, yet undiscovered, underlying principle of how nature 
organizes the neutrino sector. While the hierarchy question is of great interest, at this 
moment, there is no experimental information available. One possible exception is a hint 
from the neutrino data from supernova SN1987A [9], however the basis of this hint is under 
discussion [10]. 

To design the experiments to measure CP or T violating phase 5 we must know in advance 
#i3, or at least its order of magnitude. There have been many proposals for experiments 
which may be able to measure #13. They can be classified into two categories, long-baseline 
(LBL) accelerator experiments [11-15], and reactor experiments [16,17]. While #13 may be 
eventually measured by some of these experiments it was recognized that this measurement 
suffers from a problem of parameter degeneracy [18-23]. It stems from the fact that the 
determination of 813, 5 and the sign of Am 2 3 are inherently coupled with each other. 

A complete solution to the parameter degeneracy would require extensive measurement 
at two different energies and/or two different baselines [19,24,25], or to combine two different 
experiments [17,25-28]. We, however, seek an alternative strategy in this paper. That is, 
we propose to solve it one by one. We pursue the possibility that determination of the sign 
of Am\ 3 may be carried out with only a limited set of measurement. It is similar in spirit 
to the one employed in [21] in which we proposed a way to circumvent at least partly the 
problem of (#13, 5) degeneracy. 

In this paper, we explore the line of thought of comparing measurement at two different 
energies and/or combining two different experiments to determine the sign of Am^ and 
#13. In particular, it is of concern to us how to combine the JHF [13] and NuMI Off- 
Axis [14] experiments. Apart from possibility of earlier detection in MINOS, OPERA or 
reactor experiments, the effect of nonzero 813 is likely to be determined by observing electron 
appearance events in these LBL experiments. Hence, it is of great importance to uncover 
ways by which these two experiments can compliment one another. We discuss in this paper 
which combination of modes of operation [y or v channels) and which neutrino energies 
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divided by baseline (E/L) have better sensitivity for determining the sign of Am^ and/or 
measuring 13 . 

While our treatment applies to wider context beyond application to the particular exper- 
iments, finding a way of optimizing the JHF-NuMI interplay is, we believe, a very relevant 
and urgent question for both of these projects. It is likely that the JHF experiment in its 
phase I operates in neutrino mode [13]. Given the state of the other projects, NuMI is the 
one which is most likely to become a reality. Therefore, it is important to find the most 
profitable mode of operation of the NuMI Off-Axis project. 

We start by giving a general formalism of comparing two different measurement. It is 
to illuminate the analytic structure of the problem of two experiments comparison, and 
it would be of help in understanding the physical characteristics behind the observations 
we will provide. Instead of engaging an intensive \ 2 analysis we prefer to illuminate the 
global chart of the two experiment comparison. To carry this out we rely heavily on a 
generalized version of the CP/T bi-probability plot that was introduced and developed in 
Refs. [20,29,23]. Our discussion in this paper, while not the complete story, it is sufficient 
to enable us to understand the key physics points. The experimental issues that need to be 
addressed to complete the picture are beyond the scope of this paper. 

The reader who wants to focus on the physics conclusions, in particular on our key 
observations in JHF-NuMI comparison, can go directly to Sec. Ill, skipping Sec. II. Then, 
they might want to come back to read the Sec. II for a deeper understanding. 

II. GENERAL FORMALISM FOR COMPARISON OF TWO DIFFERENT 

MEASUREMENT 

Let us combine two measurements, either v or 9, for two different energies and/or base- 
lines. For example, one measurement could come from JHF and another measurement from 
NuMI, with both v (or v) or one v and the other v. All combinations can be treated in the 
universal formalism presented here. 

We can start with two basic equations as in [23] for small sin6* 13 : 



where 9 = sin6>i 3 , and the superscripts a and (3 label the process and the experimental 
setup. The value of the parameters X a , Y a and A a depend on the process (specified by a), 
the energy and path length of the experiment as well as sign of Amf 3 and the density of 
matter traversed by the neutrino beam. In Table I we give the values for the X's, F's and 
A's in terms of the variables X±, Y± and A 13 which will be defined below. 




(1) 
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TABLE I. The assignments for the parameters of Eq.(l) depending on the oscillation processes 



and sign of Amf ; 



The functions X±, Y± and P are given by [30] 

'23 





(2) 



Y ± = ±8c 12 s l2 c 23 s 23 [—) ^— j sin (^j sin ^— ) (3) 



with 



Aij= 2E J ' and S± = |A 13 ±aL|, (5) 

where a = y/2GpN e denotes the index of refraction in matter with Gf being the Fermi 
constant and N e a constant electron number density in the earth. Some crucial properties 
of the functions X± and Y± which allow this unified notation are summarized in Appendix 
A. 

The basic equations (1) can be solved for either sign of Amf 3 as: 



where we have ignored terms of order The right equality can be used to solve for both 
sin 5 and cos 5 from which 9 is then determined. This is a straightforwardly generalization 
of the method used in [23] to obtain, analytically the degenerate set of solutions and to 
elucidate the relationships between them. Here, we will present the result of such a general 
analysis in the Appendix A. An important point, however, is that once the solutions for S 
are known, 6 is easily obtained using left hand equality in Eq. (6). 

Experimentally the most likely channels to be realized in the near future are — > v e 
and Pfj, — > & e , therefore in the next two subsections we will specialize to these channels 
with both beams being Neutrinos (or Anti-Neutrinos) and one beam of Neutrinos and one 
beam of Anti-Neutrinos. In the following section we will consider the specific experimental 
situations presented by the JHF [13] and NuMI [14] proposals. 
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A. Both Neutrinos or Both Anti-Neutrinos 



We apply the general formalism developed in the previous section to an explicit example 
of a v-v comparison of two experiments with different energies and baseline distances. At 
the end of this subsection we will give the relationship to the v-v comparison. We will 
label the two experiments J and N in reference to JHF and NuMI. However, in spite of the 
reference to specific two projects, our general discussions are valid for any pair of the LBL 
experiments on the globe and can also be used to discuss other processes but here we will 



concentrate on v.. 



We start by introducing some convenient notation. Since we will treat Am^ > and 
AmJ 3 < simultaneously we will label the various quantities below with subscripts ± to 
indicate the sign of Am 2 3 . We define 



yJ 



2X J 



sin 



C ± = 
D ± = 

AP ± = 



± 



2X1 



cos 




si 



P J - Pi 



\ 



Yf 
2X^ 



sin 
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pN 



X^ 



'AN* 
^13 



' A N 1 
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(7) 

(8) 
(9) 
(10) 



The signs and the subscript assignment come from Table 1. 

By solving (6) with these definitions, it is easy to derive (see Appendix A), for positive 
Amjj, the allowed values of (9, S) which we assign subscripts 1 and 2; 



sin <5i 2 = 

• D + 

, 1 

COS Oi 2 = 7^— 



-S + AP + ±C + ^/D + -(AP + y 



D 



+ 



-C + AP +T S + JD + -(AP + ) 



(11) 



Notice that with this choice of signs sin 2 <5j+cos 2 Si — 1. For negative Am 2 3 , we use subscripts 
3 and 4 to distinguish them from the positive Am 2 3 solutions, 



sin 5 



3,4 



D. 



cos 5 3A = — 



-S_AP_ t C-y/D- - (AP_) 2 



-C_AP_ ± S-JD- - (AP_) 



(12) 



The relative ± signs in (11) and (12) are arbitrary at this point but are the same as those 
used in the next subsection. 

The relationships between the mixed-sign degenerate solution are given by 



cos (Si - 5 S ) 



(C + C_ + S+S- 
D,D_ 



[ AP + AP_ -JD+- APIa/PL - AP 



+ {C - S+ D ~J_ +S ) (AP-V^ - + AP +V /P- - AP 2 ) , (13) 



5 



sin (5, - 6s) = { °- S+ D ° +S ~' ( AP + AP_ - - AP^D_ - AP 2 

_ (C + C_ + S+S-) 
D,D_ 



( AP_ - AP 2 + AP+ /dT^AP 2 ) . (14) 



We note that, unlike the case of i/-P comparison [29] which will be discussed in the next 
section, there is no sensible limit A^ — > A( 3 and (aL) N — > (aL) J because in this limit no 
extra information is provided. However, one can ask the question, can we obtain some useful 
relationships between CP phases of the mixed-sign degeneracy in an approximation of small 
difference in the matter effect between JHF and NuMI? The answer is yes and we will derive 
such relations in Appendix B where we formulate such perturbative treatment. 

The solutions for 9 can easily be obtained by substituting appropriate values of sin S and 
cos S into Eq.(6). For positive Am 2 3 , the values of Q\p. are obtained by using the X 3 , Y+, A{ 3 
and X%,Y?, Af 3 in Eq.(6) and for negative Am 2 3 , the values of 6^ are obtained by using 
Xl, Yj 7 , -Af 3 andX5, , -A^. See Eq. (38) of the Appendix A for explicit expressions. 

For the v-v comparison: for positive Am 2 3 , the values of are obtained by using the 
X 3 , — Yj 7 , — A{ 3 and X^ , — Y^ , — A± 3 in Eq.(6) and for negative Am 2 3 , the values of 
#3,4 are obtained by using X 3 , — Y 3 , A{ 3 and X+ , — Y+, A± 3 . Thus the allowed region 
for the v-v comparison is identical to the v-v comparison except for the fact that the roles 
of Am 2 3 > and Am 2 3 < are interchanged. 1 



B. One Neutrinos, One Anti-Neutrinos 



We want to treat the v-v comparison in an entirely analogous fashion as the v-v com- 
parison in the previous subsection. Toward this goal we introduce a similar notation: 
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ap; 



Yl 
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XL 

2Xl 



sin 



cos 




sin 



cos 



d' ± s (s^ 2 + (c ± r, 
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(15) 

(16) 
(17) 

(18) 



where the subscript ± in (18) denote the sign of Am 2 3 and here P is for anti- neutrinos. 
The signs and the subscript assignment come from Table 1. These definitions are almost 
identical to the v-v definitions but not quite, the ± labels on the left hand side are different 
and these difference lead to markedly different outcomes. 



The difference in the sign of the Y coefficients can be compensated by taking 8 — > 5 + it, see 
Eq.(l). Thus the allowed regions are identical. 
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It is easy to verify that by using the above notation sin Si as well as 9i are given by exactly 
the same expressions as in the previous sections, e.g., (11) and (12), apart from replacement 
S± — > S'±, C± — > C'±, and so on. 

The relative ± signs in (11) and (12) are determined such that it reproduces the pair of 
the degenerate solutions, 5 3 = Si + n (mod. 2ir) in the limit of P N = P J , A J = A N and 
(aL) J = {aL) N for a given value of #13 [23]. The relation, however, receives corrections due 
to difference in matter effect between JHF and NuMI even at A J = A N . It is not difficult 
to compute the correction at P N = P J in the form 

8 3 = 5, + tt + C N ~ J [(aL) N - (aL) J \ , (19) 

to first order in JHF-NuMI difference of matter effects. The coefficient C N ~ J is given in 
Appendix B. 

Notice that the relation between the same-sign solutions of the CP violating phase is 
given generically by the formula (42) in comparison of any two channels. Among other 
things, 62 = tt — Si always holds at oscillation maximum, A N = A J = tt, despite the 
difference in matter effects at JHF and NuMI. It must be the case since the bi-probability 
trajectories shrink to straight lines. 2 



III. JHF VERSES NUMI COMPARISON 

In this section we will discuss explicitly the JHF and NuMI possibilities. JHF will have 
a baseline of 295km, and the neutrino energy at which oscillation maximum occurs is 

0.60 GeV ( ^\ 0) • (20) 

\2.5 x lO- 3 eV 2 ) K J 

For our figures we will use this energy plus 0.8 GeV which is 33% higher and near the peak 
in the event rate for v e appearance due to the rising cross section but falling oscillation 
probability assuming the same flux can be maintained. 

So far the path length for NuMI is undecided but it will most likely be between 500 and 
1000 km. We use 732km for our figures, the path length of NuMI/MINOS. At this distance 
oscillation maximum occurs at 

15 GeV ( 2.5 xffi cV j (732^) ■ < 21 > 

The other energy used is 33% higher at 2.0 GeV near the event rate peak for the reasons 
stated for JHF. 

In Fig. 1 we have plotted the allowed region in bi-probability space assume no knowledge 
of 9 and S in a comparison between JHF neutrinos and NuMI neutrinos. This allowed region 
forms two narrow "pencils" which originate at the origin and grow in width away from the 



2 The straight-line CP trajectory can be achieved even when energy distribution of neutrino flux 
times cross section is taken into account [21]. 
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origin. Associated with each of these "pencils" is the sign of AmJ 3 , thus if the pencils are 
well separated then a comparison of these two experiments can be used to determine the 
sign of Amf 3 . But the size of #i 3 is poorly determined in such a comparison. More details 
on this v-v comparison will be given in the next subsection. (A v-v comparison is identical 
but with the sign of Am^ flipped.) 

In Fig. 2 we have plotted the allowed region in bi-probability space in a comparison 
between JHF neutrinos and NuMI anti-neutrinos. Here the allowed regions are very broad 
and there is significant overlap between the allowed regions for the two signs of Am 2 ^. 
However, the size of #13 can be determined with reasonable accuracy in such a comparison. 
Fig. 3 is similar to Fig. 2 but here the concentration is on the overlap region. More details 
of this v-v will be given in the subsection following the v-v comparison. (A v-v comparison 
is identical to a v-v but with the sign of Am^ flipped.) 



A. Both Neutrinos 

It is worthwhile to have a simple formula which relates between two v^ — > v e appearance 
probabilities obtained by JHF and NuMI both in neutrino (or anti-neutrino) channel. Let 
us restrict, for definiteness, the following discussion into the case of neutrino channel. Then, 
one can easily derive the following expression 



pN _ pN 



pj p.J 

_ ^0 



xj 



cos 



sm 



'A AT 
^13 



'AW 
^13 



Yi 

xl 



cos 



sm 




= 2(C±cos5±S±sm5)6 



6 cos 8 



6 sin 5 



(22) 



We want to understand the behavior of two loci corresponding to positive and negative 
in Fig. 1 as a function of A J and A^. In particular, we are interested in how the 
slopes of the two "pencils" change in Fig. 1. To this goal we compute the slope ratio of the 
central axis of positive Amj 3 to the negative Amf 3 . Since the central axis of the "pencil" is 
obtained by equating LHS of Eq. (22) to zero, the slope of P N — P J line in Fig. 1 is simply 



given by a± = -yj-. Thus the slope ratio can be written as 



~x~l 
xJ 



sm 



B 1 ' 



sin 



x 



2 /(2 



^( B 4)/( B 4y 



(23) 



Therefore, the behavior of the slope ratio is controlled by a single function sin 2 x/x 2 . 

To have a qualitative understanding we perform perturbation expansion assuming matter 
effect is small, aL <C A13. Noting that B± = A i3 ± aL for both JHF and NuMI we obtain, 
to first order in matter effect, 



1 + 2 



2 -cotY^' 



AN 
^13 



(aL)\ 



NuMI 



- 2 



2 -cctf^ 



A/ 3 



(aL)\ 



JHF 



(24) 
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Therefore, the slope of positive-Am^ "pencil" is larger than negative- Am^ "pencil" because 
of the larger matter effect in NuMI. 

To understand the energy dependence of the slope ratio we note that | — cot a; is a 
monotonically increasing function of x in < x < ir. Let us first examine the case of 
A^ 3 = A{ 3 . When the neutrino energy is increased from the oscillation maximum (which 
means lowering A^'s) the slope ratio decreases, in agreement with the behavior in Fig. 1, 
from 1(a) to 1(d). We have checked that lowering the energy (larger A^'s) in fact leads 
to larger slope ratio. However at energies below oscillation maximum the probabilities and 
event rates rapidly get smaller. 

The behavior of slopes of "pencils" in other parts in Fig. 1 can also be understood in a 
similar manner. From Fig. 1(a) to 1(b) the slope ratio decreases, whereas from Fig. 1(a) to 
1(c) it increases. The change in 1(a) to 1(b) is obtained by lowering A^ 3 while keeping Af 3 
fixed. Since the second term in the RHS of (24) decreases as A^ decreases, and hence the 
slope ratio decreases, in agreement with Figs. 1(a) and 1(b). Similarly, from Fig. 1(a) to 
1(c) A{ 3 decreases while keeping A± 3 is kept fixed. It should lead relatively small amount 
of increase in the slope ratio, again in agreement with Figs. 1(a) and 1(c). 

However the slope of these "pencils" is not the total story, the width of the "pencils" 
is also important. Although the ratio of slopes is slightly larger in Fig. 1(c) than Fig. 1(a) 
the width of the "pencils" is significantly larger for Fig. 1(c) than Fig. 1(a) such that the 
separation of the allowed regions is smaller for Fig. 1(c) than Fig. 1(a). The square of the 
width of the "pencils" is controlled by the quantity 

MSl + CD (f ) (|) cos(f - M). (25) 

For A^ = A/ 3 , the width equals 




X. 



N 




(26) 



and is very small. This smallness follows from the fact that, at the same A 13 , the identity 

V N V J 



= holds ( depends only on vacuum parameters), and ratio of X± to X± is 



Y N Y J 

close to unit. If A^ 3 ^ A J l3) then —7= 7^ 7= and the width of the pencil grows rapidly 

V- v :. v x ± 

as the cancellation that occurs for the same A13 no longer holds. Thus the width of the 
"pencils" in Fig. 1(b) and Fig. 1(c) are approximately equal and are much larger than the 
width in Fig. 1(a) and Fig. 1(d). 

The conclusion to be gained from these results and figures is that if you are making a 
comparison of a JHF neutrino experiment and a NuMI neutrino experiment to determine 
the sign Am\ 3 then the best separation occurs when A J = A^ i.e. the same || for both 
experiments. Smaller values of the E are slightly preferred so that we expect the optimum 
value, once all experimental issues are included, to be near oscillation maximum (Fig. 1(a) 
and 1(d)). Away from the same ^ for both experiments the separation between the allowed 
region becomes worse, apart from special ranges of the CP or T violating phase 5 when 
(f ) \jhf > (f ) \nuMI (Fig. 1(c)). Whereas for the choice (f ) \ JHF < (J ) \ Nu mi (i-e 
A J > A^) there is significant overlap between the two different sign of Am\ 3 allowed regions 



9 



(Fig. 1(b)). (If one moves further, say, to E JHF = 0.6 GeV and E NuM i = 2.5 GeV, the two 
pencils overlap almost completely.) In these regions no matter how accurate the experiment 
the sign of AmJ 3 cannot be determined by this comparison. (For JHF anti-neutrinos versus 
NuMI anti-neutrinos the conclusion is the same.) 

In summary, for the comparison JHF neutrinos to NuMI neutrinos to be useful in the 
determination of the sign of Amf 3 , choosing 

JHF > ( -r) \NuMI (27) 



L 

is of great importance with the preference for equal (j^j- For \jhf < ) \nuMI there 
is significant overlap between the narrow allowed regions so that this is not a good choice 
for this comparison. 

Contrary to advantage of v-v comparison in determination of the sign of Am^, mea- 
surement of #13 suffers from large uncertainty in this channel. The bi-probalitity trajectory, 
for a given 9, moves nearly along the direction of the "pencil" as one varies 5, thus any 
measurement of the oscillation probabilities with finite resolution cannot pinpoint the value 
of 9. To indicate this point we have plotted in Fig. 1 by numbers in % beside thin arcs 
the fractional difference A9/9 = (9 2 — 9\)/ ((9 2 + 0i)/2), where 9 2 and 9\ here indicate the 
maximal and the minimal values of #13 along each arc. The numbers are large, typically 30 
% or even larger. 

In fact, it is easy to do analytical estimate of the fractional difference. Suppose that 
we have measured an appearance probability P J or P N . We assume for simplicity that 
we know the sign of Am^, and restrict ourselves into the case of oscillation maximum, 
A{ 3 = A^ 3 = 7r. The measured probability allows a range of #13 [21] and one can easily 
calculate the maximal value of A9/9 within the range, which leads to 

A9 n A12 /sin(aL/2)\ , . 

Since P J / N is a few % level, it leads to a few tens in % of the fractional difference, in 
agreement with the numbers in Fig. 1. 



B. One Neutrinos and One Anti-Neutrinos 



We now turn to the comparison for JHF neutrinos and NuMI anti-neutrinos or vice 
versa with flipped signs of Amf 3 . Again one can write a simple formulae relating the two 
oscillation probabilities similar to Eq. (22) 



J AT 



pN 



pj pj 



VJV 



COS 



'AW 
^13 



xi 



cos 



sin 



^13 



xi 



sin 




= 2( y C±cos5±S±sm5)9 



9 cos 5 



6>sin<5 



(29) 



The square of the RHS of this equation controls the square of the width of the allowed region 
and is given by 
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4(S| + C«)= ^ + ^ + 2 ^ P s (? + ?) (30) 



apart from an overall factor of 9 2 . As (A^+A^) varies from n (both experiments at energies 
twice the oscillation maximum energy) to 2n (both experiments at the oscillation maximum 
energy) this squared width grows from 






to - • (31) 



Note that the later width is larger than the first width since Y_/X_ has opposite sign to 
Y + /X + (no cancellation occurs here). Thus the width of the allowed region grows as you 
lower the energies of the experiments to oscillation maximum. This can be seen in Fig. 2; 
2(d) has the smallest width whereas 2(a) has the largest in accordance with the above 
statement. 

There is significant overlap between the allowed regions of positive and negative Am^ 
contours. Hence, the v-v comparison does not appear to be the right way for determining 
the sign, unless 5 turns out to be close to 37r/2 and it/ 2 for Amf 3 > and AmJ 3 < cases, 
respectively [20]. However, this comparison may be better to measure 9 13 if JHF operate only 
in neutrino channel, as planned, and if NuMI Off-Axis is running with significant overlap 
with JHF phase I. For this purpose, the best way of measuring # 13 would be to tune the 
energy at oscillation maxima for both JHF and NuMI. This is just an extension of the KMN 
strategy [21]. 

The reader may be curious as to why the slopes of the shrunken trajectories at A 13 = n 
are positive in v-v comparison whereas they are negative in v-v comparison. This can be 
easily understood as follows: The slopes of shrunken trajectories in v-v comparison are 
given by ' /Y+ (Y N /Y^) for positive (negative) Am^. On the other hand, the slopes 
in v-v comparison are given by Y^ jY+ {Y+ /Y^) for positive (negative) Amf 3 . Since Y± 
differ in sign, as in Eq. (3), the latter slopes are negative definite, whereas the former are 
positive definite. The fact that the v-v trajectory is along the "pencil" and v-v trajectory 
is perpendicular to axis of the "cigar" can be understood by a similar argument. 



IV. ANOTHER INTERESTING COMPARISON 

When we have access to v e and v e beams then a comparison of a v^ — ■> v e experiment 
with v e Vfj, experiment will be possible. 3 This comparison is interesting because it directly 
compares two CPT conjugate processes. If CPT is conserved, then at the same E/L the 
only difference between the oscillation probabilities for v^ v e and v e — > v^ can come 
from matter effects. Thus this comparison is even more sensitive to the mass hierarchy, i.e. 
the sign of Am^, than the neutrino- neutrino comparison discussed earlier in the previous 
section. 



3 Possible ways to produce electron neutrino beams include muon storage rings (Neutrino Factories 
[31]) and /3-beams [32]. 
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For this comparison the bi-probability figure looks similar to Fig. 1 but the difference 
in slope between the "pencils" is enhanced because the sign of the matter effect is opposite 
for V e — > Dp compared to — > v e and the CP or T violating terms are identical in sign and 
magnitude. Therefore this comparison provides an excellent way to separate the two signs 
of Amy. If we label the two experiments as J and N as before then the ratio of slopes is 
given by 



xl 



AN 
^13 



cot 



AT 
43 



(aL)\ N + 2 



A/ 3 



cot 



'A 



13 



(32) 



This expression differs from that of Eq.(24) by the sign of the third term. Here the matter 
effects for both experiments enhance the ratio of the slopes whereas in the neutrino-neutrino 
comparison early there was a partial cancellation. Although at present the comparison of 
this section is purely academic it is instructive and useful for understanding the general 
nature of the comparison of neutrino oscillation experiments and will be important to test 
the possibility of CPT violation in the future. 



V. SUMMARY AND CONCLUSIONS 

We have presented a general formalism for comparing two or more long baseline neutrino 
oscillation experiments and applied this formalism to the JHF and NuMI experiments. The 
combination of modes that will be important depends on the question one is asking and 
what other information is available at the time. The use of bi-probability diagrams like 
the ones presented in this paper will be important for understanding the physics issues and 
making trade offs at the time decisions are made regarding neutrino energies and modes of 
operation. 

In general, the both neutrino or both anti-neutrino comparison is useful for determining 
whether the mass hierarchy is normal or inverted, i.e. the sign Amj 3 . However, it is impor- 
tant here that the experiment with the larger matter effect (longer baseline - NuMI) have 
a smaller or equal neutrino energy over baseline, E/L, than the experiment with smaller 
matter effect (shorter baseline - JHF). If the JHF experiment runs first then NuMI should 
also run at a similar or smaller E/L to provide the best sensitivity to the different mass hi- 
erarchies. The separation between the different mass hierarchies is reduced if the NuMI E/L 
is higher than the JHF E/L. The size of the small mixing angle, #13, cannot be determined 
to better than about 30% with this comparison alone. 

The one neutrino - one anti-neutrino comparison is in general most useful for determining 
the size of the small mixing angle, 6>i 3 and the CP or T violating phase 8. It could also be 
useful in determining the mass hierarchy provided nature does not choose the large overlap 
region in the bi-probability plot. For this comparison, the uncertainty in #13 from the 
degeneracy issue is not important until the experimental resolution in 9 13 is better than 
10%. 



12 



ACKNOWLEDGMENTS 

HM and HN thank the Theoretical Physics Department of Fermilab for warm hospitality 
extended to them during their visits. This work was supported by the Grant-in-Aid for 
Scientific Research in Priority Areas No. 12047222, Japan Ministry of Education, Culture, 
Sports, Science, and Technology. Fermilab is operated by URA under DOE contract No. DE- 
AC02-76CH03000. 



APPENDIX A 

In this Appendix, we first explain briefly why the universal formalism with generic nota- 
tions given in Table 1 is possible. The first secret behind it is the relationship between X's 
and Y's with different sign of Am^. Let us think of the expressions of the appearance prob- 
abilities in neutrino and anti-neutrino channel similar to (1), and define coefficient function 
X± as well as X±. They can be expressed by using a single function X as X± = X(±Am 2 , a) 
and X± = X(±Am 2 , —a). The similar notation can also be defined for Y's. Then, 



X ± =X T 
Y± = —Y T , 



(33) 



because X's and F's (except for extra sign) are the function only of A 13 ± aL, which follows 
from the CP-CP relation in [29] and the approximation introduced in [30]. We note that 
X± and Y± satisfy the useful identity 



+ 



TxT' 



(34) 



Next we present general solutions of (1) whose validity extends any combinations of 
channels tabulated in Table 1. For any given process and sign of Am^ there are two 
solutions to these equations which we will label by (9\,Si) and (02, 62)- The equations (6) 
can be solved for the CP phase 5 for either sign of Am 2 3 as: 



sm#i ; 2 = 
cos Si 2 = 



-SAP ± C^S 2 + C 2 - (AP) 2 
-CAP t StJs 2 + C 2 - (AP) 2 



/(S 2 + C 2 ) 
/(S 2 + C 2 ) 



(35) 



Notice that with this choice of signs sin 2 5i + cos 2 5j = 1, and the variables AP, S, and C 
are defined as 




S = 



C = 



A a \ (A 13 * 
/?° sin ( — J /,'Sin I — 



-R a cos (—) + R P 

cos — 



(36) 
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where we have defined the variable 



R a = 7^ (37) 
2X a y ' 



which appears frequently. The corresponding 6>'s are then given by 



A a a A 13 ( A a A f3 \ \ 
+ R a sin {~2~)R P sin (— ) I cot — - cot — J cos 6 1>2 J 

/ A a A 13 \ 

/(P*sin( — )-P^sin( — )J. (38) 

Under the interchange of a and (3, Si <-> S 2 and 9\ <-> #2 as it should. Also for useful exper- 
iments cot 4r is always finite. These two solutions coincide, 61 = 9 2 , when the coefficient 
in front of the cos5i j2 term vanishes or cos^i = cos5 2 - The first possibility occurs when 
A a = A 13 or A a = ±7r = ±AA While the second occurs at the edge of the allowed region 
when the square root in Eq.(35) vanishes. 

The allowed region in (PA P a ) bi-probability space is given by 

/A" A^\ 

(AP) 2 < S 2 + C 2 = {R a f + (P' 3 ) 2 - 2R a R 13 cos I J . (39) 

This region is determined by requiring the sin<5 and cos 5 to be real. The boundary is 
determined by the equality in eqn (39). Thus the quantity 

/A« A/3\ 

{R a f + (P^) 2 - 2R a R? cos I =- - =- J . (40) 

controls the separation of the two boundaries and hence the size of the allowed region. 
Depending on the value of (A a — A' 3 ) this can range from (AP) 2 = (R a — R 13 ) 2 when 
(A a -A 13 ) = 0, 2?r • • • to (AP) 2 = (R a + R 13 ) 2 for (A a - A' 3 ) = ±tt, ±3tt • • •. Which of these 
two extremes gives the larger allowed region depends on the relative signs of R a and PA 
This will be very important since the relative signs changes as we go from both neutrino (or 
both anti-neutrino) to one neutrino and one anti-neutrino process. 

Let us now compute the difference between two solutions of S. One can show by using 
(35) that 



P«sin(4r) -P^sinf^ 

cos (<$i + <S 2 ) = 1 - 2 l - ^ y A /J — T . (41) 

1 ' (P-) 2 + {R?y - 2R«R? cos - f ) K ' 

One can easily check that Eq. (41) reduces to Eq. (45) of [23] if we take A a = -AA The 
relationship (41) implies that 



\R a sin sinf^)| 
5 2 = TT - 5i + 2 arcsin V 2 ' V 2 ' (42) 



(P«) 2 + {RI 3 ) 2 - 2R a R' 3 cos - f- 
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where use has been made of the relation arccos (1 — 2a; 2 ) = 2arcsina;. Thus the second 
solution 82 differs from 7r — 81 by a constant which depends on the energy and path length 
of the neutrino beams as well as on the process but not on the mixing angle #13 for any 
two-experiment comparison. This generalizes the result obtained in [23]. 

At the oscillation maximum, |A Q | = {A 13 ] = n and cos (8\ + 8 2 ) = —1 i.e 8 2 = n — 81. This 
is a reflection on the fact that at oscillation maximum the trajectory in the bi-probability 
space is an ellipse with zero width as can easily be seen from Eq. (1) as there is no dependence 
on cos 8 at this point. 



APPENDIX B 

We evaluate the correction to the relationship between two degenerate solutions of 8 in 
v-v and v-v comparison. Let us first discuss v-v comparison first. 

We expand S' ± , C' ± , and AP± at around the JHF parameters. For simplicity, we only 
deal with the case A N = A J . They read, 



± 



4 

e_A (Y 
4 



n, 



cos 



' A J s 

^13 



9 



' aL 



2 j 

'aL r 



T9 
T9 



API = AP£ t — 



2\ 



P N -Pg (B± 

xi 9 \Y 



where €a = (aL) N — (aL) J . We have defined the function g as 



B± 
2 

B± 
2 



(43) 
(44) 

(45) 



q(x) = cot x. 

v 1 x 



(46) 



The function g monotonically increases in < x < ir. 

Then, sin (81 — 8 3 ) in Eq. (14) can be expanded to first order in as sin (81 — <5 3 ) 
C N ~ J 6a- The coefficient C N ~ J is given by 



C N - J - S t CJ g( aLJ ) I 1 AP ° 



+ 



D J 
1 



■ 2 J ' 4 - APi : 

S J C J AP n 



9 



2sinA^ D J ^J D J _ AP 2 [ 



9 



+ 9 

'B J ' 



'B^ 



^9 



aL" 



1 sin A" 
8 D J 



'XT 



B J _ 



+ 



2 y'/V - API 



P N _pN 



(B J \ 



'Bi 



P N _pN 



(47) 



where S J , C J , and AP indicate, respectively, the 6a — > limit of S'±, C'±, and AP' + (or 
superscript N replaced by J). Notice that AP' + = —AP'_ = AP in the limit P N = P J and 
that S' = Si = S J andCi = Cl = C J . 
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Now we turn to v-v comparison. In this case there is no zeroth order term in S± etc. for 
reasons explained in subsection IIA. Instead they start with the first-order terms in e^. We 
do not give the details but just mention that sin (Si — S 3 ) is given to first-order by 



sin (5 1 - S 3 ) 



x 



y-JyJ 

V x + X - 



sin ( A J 



2 



1 + 



(y/) 2 

4XiP J 



9 



aL J 
2 



B J _ 

9\ — 



cos (A J 



{Y J } 2 



AXiP J 



i+ Ml 



B J _ 

9\ — 



\ 



(Y- 



J\2 



AXlP J 



9 



aJJ_ 
2 



9\~f 



- 1 - 



(Y- 



J\2 



AX J _P J 



9 



2 



9\~f 



- 1 



One can show that (Si — S 3 ) is in the second quadrant because cos (Si — S 3 ) < in our 
convention. 
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FIGURES 



A0/0 (%) for positive and negative Am 13 
JHF neutrino vs. NuMI neutrino 
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FIG. 1. Allowed range of P{yu — > v e ) for JHF verses P{y^ — > v e ) for NuMI, which are referred 
to as "pencils" in the text, are delimited by thick solid (dashed) lines for positive (negative) 
Am 2 3 for the energies (Ejup/GeV, Ej^ u mi /GeV) = (a) (0.6,1.5), (oscillation maximum for both 
experiments) (b) (0.6,2.0), (c) (0.8,1.5) and (d) (0.8,2.0). In the same plot, the positions for some 
representative values of the fractional variation across the width of the "pencil" of 6 = sin #13, A6/6 
[%], indicated by numbers, are shown by thin solid arcs. Inside each allowed region, trajectories 
corresponding to sin 2 2#i3 = 0.02, 0.05 and 0.09 are plotted by dotted lines. The mixing parameters 
are fixed to be |Amf 3 | = 2.5xl0~ 3 eV 2 , sin 2 26» 23 = 1.0, Amf 2 = +7xl0~ 5 eV 2 and sin 2 26 12 = 0.85 
whereas #13 and 5 are assumed to be unknown. The electron density is fixed to be Y e p = 1.15 and 
1.4g cnr 3 for JHF and NuMI experiment, respectively. For JHF and NuMI both anti-neutrinos, 
the roles of Amf 3 > and Am 2 3 < are interchanged. 
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A0/0 [%] for positive and negative Am 
JHF neutrino vs. NuMI anti-neutrino 
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FIG. 2. Iso-contours of A8/9 [%] are indicated in the allowed range of -PA M — * t> e ) for JHF 
verses P {pa — > z? e ) for NuMI delimited by solid (dashed) lines for positive (negative) Am? for the 
same combinations of energies as in Fig. 1. Inside each allowed region, trajectories corresponding 
to sin 2 2#i3 = 0.02, 0.05 and 0.09 are plotted by dotted lines. The assumption about the mixing 
parameters are the same as in Fig. 1. For JHF anti-neutrinos and NuMI neutrinos, the roles of 
Amf 3 > and Amf 3 < are interchanged. 
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A0/0 [%] in the overlap region of different signs of Am 
JHF neutrino vs. NuMI anti-neutrino 
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FIG. 3. Iso-contours of A9/9 [%], except for (a), are indicated in the overlap regions of the 
allowed range of P(va — > u e ) for JHF verses P (pn — > v e ) for NuMI for the positive and negative 
signs of Amfg for the combinations of energies shown in Fig. 2. For (a), inside the shaded region, 
A9/6 varies very little, taking values close to 9 %. 
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